Spectral analysis of topological defects in an artificial spin-ice lattice 
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Arrays of suitably patterned and arranged magnetic elements may display artificial spin-ice struc- 
tures with topological defects in the magnetization, such as Dirac monopoles and Dirac strings. 
It is known that these defects strongly influence the quasi-static and equilibrium behavior of the 
spin-ice lattice. Here we study the eigenmode dynamics of such defects in a square lattice consisting 
of stadium-like thin film elements using micromagnetic simulations. We find that the topological 
defects display distinct signatures in the mode spectrum, providing a means to qualitatively and 
quantitatively analyze monopoles and strings which can be measured experimentally. 



Spin ices are magnetic structures in which the inter- 
actions lead to frustration resulting in complex mag- 
netic ordering and collective behavior. Examples of 
such spin ice systems on the atomic scale are the py- 
rochlore structures [T]-[3] Dy2Ti2 07 and Ek^T^Oy. In 
these, the atomic spins arranged at the vertices of tetra- 
hedra are forced by crystal fields to point either in or 
out of the tetrahedra. Artificial spin-ice lattices are 
patterned nanoscale structures in which each element 
is small enough to have essentially uniform magnetiza- 
tion. The geometrical arrangement of the structures in 
the lattice leads to frustration by design. Examples of 
artificial spin ice are the Kagome lattice [4j |5] and the 
stadium square lattice, consisting of stadium-shaped el- 
ements on a square lattice [SHE]- The magnetization at 
the vertices of the stadium square lattice satisfies simple 
spin-ice rules shown in Fig. [I] The equilibrium configura- 
tion and magnetization reversal of such lattices have been 
the subject of intense study [9-12 . Topological defects, 
in which the magnetization at a vertex is in a one-in- 
three-out (or three-in-one-out) configuration have been 
observed[9l [TTJ [13-15 . In such cases, the flux integral 
over a surface enclosing the inner ends of the elements 
at a defect vertex is non-zero, indicating that the ver- 
tex contains a net magnetic charge, similar to a Dirac 




Figure 2. (Color online) a) Simulated singly-charged 
monopole pair (G + G~) linked by a Dirac string (DS). Vor- 
tices V form along the string. The reference state of the lat- 
tice is shown in pale blue, while the string and monopoles are 
highlighted, b) Each vertex with a doubly- charged monopole 
(G* + or G*~) requires the presence of an oppositely-charged 
monopole (G~ or G + ). The colormap corresponds to the 
projection of the magnetization along the (11) direction of 
the string. 
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Figure 1. (Color online) Different vertex configurations in a 
square spin ice lattice. The first two (RS and DS) are charge 
neutral while G and G* show two possible charged configu- 
rations. The arrows and colors represent the magnetization 
direction and the orientation with respect to the vertex: or- 
ange "m" and blue "out" . The magnetic charge is defined as 
positive when the net magnetization points out of the vertex. 



monopole [16 (see Fig. [T]). Such monopoles typically oc- 
cur in pairs (e.g. monopole- antimonopole) connected by 
a defect string, or "Dirac string" along which the mag- 
netization is reversed with respect to a lattice reference 
state (Fig. |2|. Recent studies have shown that these de- 
fects profoundly affect the equilibrium behavior and the 
magnetization reversal of such lattices [5j [9j [TTJ [12j [T7l 
[19] • An interesting, yet unexplored question is whether 
the presence of such defects affects the resonant dynamics 
of the spin-ice lattice. Here, we investigate the frequency 
spectrum of eigenmodes, or resonances, of topological de- 
fects in a square spin ice lattice using micromagnetic sim- 
ulations. We show that each topological defect displays 
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Figure 3. (Color online) (a) Simulated lattice composed of 112 magnetic elements. It contains four monopole-antimonopole pairs 
connected by Dirac strings extending over 28 elements, (b) Evolution of the magnetization dynamics spectrum with increasing 
string length and number of monopole-antimonopole pairs compared to the reference state. The shaded labels correspond to 
the main, distinct signatures of the topological defects. 



distinct and localized features, both spatially as well as in 
frequency. These features act as a fingerprint for the type 
of defect and moreover allow to quantify the number of 
defects present in the lattice. We also show how spectral 
features in a lattice arise from the magnetostatic coupling 
between the edge-like modes of isolated elements. The 
coupling leads to multiplets that are split depending on 
the local magnetization configuration. In addition, the 
fundamental uniform mode of a single stadium is split in 
the presence of topological defects. 

We have studied the resonances of square spin ice 
lattices consisting of 112 Permalloy (Py) elements 
(Fig. [3^t) with micromagnetic simulations using a fully 
three-dimensional finite-element algorithm based on the 
Landau-Lifshitz-Gilbert equation [20 . Each element 
(stadium) is 290 nm long, 130 nm wide, 20 nm thick 
and is discretized into about 22800 tetrahedral elements, 
corresponding to a cell size of about 3 nm 3 . The lattice 
constant of the array is 390 nm. The material parame- 
ters chosen for Py are the following: A = 13 pJ/m (ex- 
change constant), [ioM s = 1.0 T (M s : saturation mag- 
netization), and K u = (magnetocrystalline anisotropy) 
[2T] . We prepared a reference state (RS) of the lattice 
by saturating the array along the (11) direction and adi- 
abatically removing the external field. The system then 
relaxes into a state with the magnetization oriented along 
the +y direction in the vertical elements and along the 
-fx direction for the horizontal elements of the array. We 
insert topological defects with charges ±1 and ±2 (de- 
noted and G* ± , respectively) [22] connected by Dirac 
strings (denoted DS) of variable length by reversing the 
magnetization in selected stadiums, and letting the mag- 
netization relax to an equilibrium state (see Fig. [2|. The 
dynamics and spectrum of eigenmodes of the lattice is 
then obtained by applying a short perturbation (10 mT 
field pulse in the (11) direction for 20 ps) and integrating 
the LLG equation in time. 

The excitation spectrum of the RS is shown in Fig. [3^b) 
as the grey, filled area. It is characterized by a few promi- 



nent peaks corresponding to distinct eigenmodes. The 
first one (at position (T)), slightly above 2 GHz, corre- 
sponds to edge excitations of the outermost stadium ele- 
ments (rim) of the lattice. A weaker mode appears above 
4 GHz corresponding to edge excitations of the mag- 
net ost at ically coupled elements inside the lattice. Two 
modes, labeled (4), are visible slightly above 5 GHz, in 
which the rim elements display higher order edges modes 
coupling to the inner elements of the lattice near the rim 
elements. Above 6 GHz, label (F) indicates a mode corre- 
sponding to higher-order edge oscillations of the individ- 
ual elements. This mode is split: owing to their smaller 
magnetostatic coupling, the lattice rim elements oscillate 
at a slightly lower frequency than the inner elements. Fi- 
nally, a mode at 8 GHz (position (6)) corresponds to the 
ferromagnetic resonance in which the oscillation is ap- 
proximately uniform in the interior of the elements. The 
colored lines in Fig. [3^b) depict the spectrum evolution 
as topological defects and Dirac strings are inserted into 
the lattice. The RS spectrum is clearly modified with 
some mode frequencies shifting, and new peaks appear- 
ing. Near mode (T) of the RS a new mode associated with 
oscillations of the G + G~ pairs appears. The amplitude 
of this mode increases with the number of monopole- 
antimonopole pairs. Later, we show that this mode is 
close in frequency to the edge modes of the RS, but is 
spatially localized at the monopoles. Therefore, it can be 
distinguished from the edge modes either in the thermo- 
dynamic limit or by using detection techniques sensitive 
to the interior of the lattice [23 . Similarly, around 5 GHz, 
a new mode associated with doubly-charged monopoles 
unfolds alongside the RS mode at (7). Above 6 GHz, 
oscillations of the monopole-antimonopole pairs give rise 
to a new mode within the lattice (red curve in Fig. [3] at 
label (J))- When a large number of pairs is present, this 

mode becomes dominant (green curve). Finally, at 
the ferromagnetic resonance peak is split and shifted in 
the presence of Dirac strings. In addition to the mod- 
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ifications in the GS spectrum, two new resonances are 
observed: mode (2) at 3 GHz corresponds to oscillations 

of the Dirac string edges while mode (3) at 3.5 GHz is 
due to oscillations of doubly-charged monopoles. 
We now investigate in detail the modes associated with 
the different topological defects. Fig. [4^i refers to the 
mode associated with a G + G~ pair at position (T) in 
Fig. [3] The static magnetization configuration is indi- 
cated by the arrows in the first column. The color maps 
(columns two and three) depict the oscillations of the 
x and y components of the magnetization, respectively, 
which are clearly localized at the monopoles. A split- 
ting in the frequency peak is associated with the G + G~ 
resonance (graph, fourth column): Monopoles located 
at the edges of the lattice (green curves) contribute to 
the edge oscillations, thereby increasing the amplitude 
of the RS mode (indicated in the figure) slightly above 
2.2 GHz. Monopole- type configurations away from the 
edges oscillate at a lower frequency, around 2.1 GHz in 
this case (maroon curves). The peak amplitude increases 
linearly with the number of monopole pairs in both cases 
as shown in the last column (simulations were carried-out 
for a maximum number of five pairs). It is therefore pos- 
sible, not only to identify the presence of monopoles in 
the lattice, but equally to quantify the number of pairs, 
given a reference point. Fig. [4)3 refers to the oscillations 
of the stadiums in the Dirac string connecting a G + G~ 
pair position (2) in Fig. [3] The spectral amplitude of 
this peak increases linearly with the total number of ver- 
tices in the string (string length). Fig. Up (position (3) 
in Fig. |3| shows the oscillations of a aoubly-charged 
monopole pair. Again, we note the linear relationship 
between the number of monopole pairs and the spec- 
tral amplitude of the mode. From this linear relation- 
ship, it can be concluded that the different strings and 
monopoles are decoupled from each other, even though 
they are rather closely packed (see Fig. [3^1) [24 . Finally, 
Fig. 4H (position ((T) in Fig. 3 illustrates the splitting of 
the ferromagnetic resonance (FMR) mode [25] . The spa- 
tial oscillations are shown here by projecting the magne- 
tization along the (11) direction of the string, as in Fig.|2j 
in order to optimize the contrast. The two color maps 
correspond to spatial oscillations at the bulk FMR fre- 
quency of 8 GHz (second column) and at 8.3 GHz (third 
column) , showing that the splitting stems from the Dirac 
string, which oscillates at a higher frequency. The peak 
amplitude of the mode at 8.3 GHz therefore increases 
with increasing string length at the expense of the FMR 
peak at 8 GHz, as shown in the amplitude plot (rightmost 
panel). Again, a slight splitting occurs when the strings 
terminate at the lattice edges (green curves) or within 
the lattice (maroon curves) because of the difference in 
magnetostatic coupling at the string extremities. 
The magnetostatic coupling at the vertices drives the 
spectral features of the monopoles and the Dirac strings, 
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Figure 4. (Color online) Examples of local oscillations of the 
magnetization (Am) characteristic for monopoles (a), Dirac 
strings (b) and doubly- charged monopoles (c). The first col- 
umn shows the static magnetic configuration. The second and 
third columns respectively show the oscillations of the x and 
y components of the magnetization. The plots quantify the 
increase in amplitude of the spectral signature of the differ- 
ent defects after subtraction of the reference state. The last 
row (d) illustrates the splitting of the ferromagnetic resonance 
peak in the presence of a Dirac string. Labels are from Fig. [3] 



and the magnitude of this coupling depends on the lat- 
tice constant. The experimental observability of these 
features clearly depends on their robustness with respect 
to variations in the lattice constant. Figure [5] shows sim- 
ulations with different lattice constants for a given con- 
figuration composed of one monopole-antimonopole pair 



linked by a string spanning six stadiums. For compari- 
son, the spectrum of a single stadium with initial equilib- 
rium magnetization along the long axis of the element is 
shown (full shaded graph). There are three main peaks 
below 10 GHz for the single stadium. The first peak at 
1.4 GHz corresponds to two modes of even and odd sym- 
metry in which the dynamical part of the magnetization 
is localized near the edges. Similarly, the peak at 5.7 GHz 
corresponds to two degenerate higher-order edge modes 
(one even and one odd) with an additional nodal line 
near the edge. (The excitation pulse used here couples 
only to the even modes.) Linear combinations of these 
edge modes couple magnetostatically in a lattice and give 
rise to multiplets, the splitting of which depends on the 
equilibrium structure of the magnetization at each vertex 
and on the lattice constant. Some of these linear combi- 
nations are localized near the vertices and are sensitive 
to the local magnetization, and therefore to topological 
defects near the vertices. Finally, the peak at 8 GHz 
corresponds to the (bulk-like) uniform ferromagnetic res- 
onance mode. In the presence of topological defects, with 
low magnetostatic coupling (lattice constant of 450 nm) 
the original peak at 1.4 GHz vanishes as the frequency 
of the edge oscillations is blue-shifted to 2.8 GHz (in- 
dicated by the square frame in the figure). Given that 
monopole modes are connected to oscillations of the ele- 
ment edges, a peak associated to the monopole pair also 
appears at 3 GHz. At lower frequencies (1 GHz), another 
distinct peak corresponds to oscillations of the monopole 
pair. This peak is blue-shifted when the lattice constant 
is decreased, and its frequency becomes comparable to 
the frequency of the lattice rim excitation at a spacing of 
390 nm, as described earlier (Fig. 3 o)). An equally dis- 
tinct peak corresponding to the Dirac string oscillation 
appears at 2 GHz for a 450 nm lattice constant and is red 
shifted. The higher order edge mode found in a single el- 
ement (at 5.7 GHz) is robust: it persists for every lattice 
separation, eventually splitting as described earlier. As 
expected, the frequency of the ferromagnetic resonance 
mode remains unchanged. It is interesting to note that 
even for a lattice constant of 450 nm, the elements within 
a lattice remain sufficiently coupled to display features 
essentially different from the single stadium spectrum. 
In summary, we have demonstrated using micromagnetic 
modeling that topological defects in an artificial spin-ice 
lattice gives rise to distinct features in the spectrum of 
eigenmodes. These features are due to modes localized at 
the topological defects. The presence and nature of topo- 
logical defects display a one-to-one correspondence with 
the spectral features, which thus provide a means of ex- 
perimentally identifying the presence of the different de- 
fects. Moreover, the spectral amplitude of these features 
is proportional to the number of defects present, enabling 
a quantitative analysis based on the spectra. The appear- 
ance of localized modes at topological defects enables 
new possibilities of manipulating magnetization modes 
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Figure 5. (Color online) Evolution of the spectral features 
with decreasing element separation for an array containing 
one monopole-antimonopole pair connected by a string of ten 
elements, compared to a single stadium. As the lattice con- 
stant of the array is reduced, the spectral features are red- 
shifted while the amplitude of the ferromagnetic resonance 
peak decreases. G: peak characteristic of the monopole pair 
oscillations. DS: peak characteristic for the Dirac string os- 
cillations. 



in spin-ice lattices, for example by inserting monopole 
pairs using spin polarized currents and propagating Dirac 
strings through the lattice. This may allow emerging ap- 
plications combining magnonics, in which propagating 
spin waves in patterned magnetic thin films are manip- 
ulated, with artificial spin-ice lattices, in which localized 
spectral features based on topological defects can be tai- 
lored. 
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